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We show how to use dimensional regularization to determine, within the Arnowitt-Deser-Misner 
canonical formalism, the reduced Hamiltonian describing the dynamics of two gravitationally in- 
teracting point masses. Implementing, at the third post-Newtonian (3PN) accuracy, our procedure 
we find that dimensional continuation yields a finite, unambiguous (no pole part) 3PN Hamiltonian 
which uniquely determines the heretofore ambiguous "static" parameter: namely, u s = 0. Our 
work also provides a remarkable check of the perturbative consistency (compatibility with gauge 
symmetry) of dimensional continuation through a direct calculation of the "kinetic" parameter u>k, 
giving the unique answer compatible with global Poincare invariance (wj. = ||) by summing ~ 50 

■ different dimensionally continued contributions. 

(N 

I. INTRODUCTION 

The problem of the gravitational interaction, within Einstein's theory, of two bodies has been revived by the 
realization that the most promising candidate sources for ground-based interferometric gravitational-wave detectors 
such as LIGO and VIRGO are binary black holes. The most significant part of the gravitational-wave signal emitted 
,_H , by two black holes corresponds to the last few orbits before coalescence. This makes it urgent to derive the equations 

■ of motion of two black holes with the highest possible accuracy. The equations of motion of two compact bodies (black 
holes or neutron stars) at the 2.5 post-Newtonian (2.5PN) approximation were first obtained in Refs. [We recall 

• that "nPN approximation" refers to the terms of fractional order (v/c) 2n ~ (GM/(c 2 r)) n in the equations of motion.] 
Of special importance for the following was the proof || that the property of "effacement" of the internal structures 
of (non-spinning) gravitationally interacting compact bodies implied that their equations of motions depended only 
| on their two "Schwarzschild" masses mi, up to the 5PN level. This effacement property was used to prove in (J| 
that the 2.5PN equations of motion of two compact objects were correctly obtained by: (i) modelling each compact 
""^ " object by a point-mass, i.e. a Dirac-delta-function source, and (ii) regularizing the divergencies generated by the use 
qi of delta- function sources by Riesz' analytic continuation method [Q. It was also mentioned in || that the correct, 
"effaced" equations of motion could as well be obtained by regularizing the 2.5PN divergencies by using dimensional 
regularization. 

More recently the 3PN equations of motion of two point-masses have been obtained by two separate groups: 
• IpH ! lH~U anc ^ 110' usm S different approaches, different coordinate gauges, and different regularization methods. The 
first group used the Arnowitt-Deser-Misner (ADM) canonical approach jy], and derived the two-body Hamiltonian 
corresponding to the ADM transverse-traceless (ADMTT) gauge. It regularized contact terms by Hadamard's partie 
finie, and divergent integrals by a combination of a Riesz-implemented Hadamard's partie finie approach with (when 
applicable) ordinary distribution theory |p|,^2|J^] . The second group worked with the harmonically relaxed Einstein 
equations and derived the equations of motion in a harmonic gauge. It regularized the contact terms and the divergent 
integrals by means of special-purpose variants of Hadamard's partie finie approach applicable to a certain class of 
"pseudo functions" (f^Jl^]. In spite of the considerable efforts spent by both groups, their regularization methods 
are defective in two (interconnected) respects: (i) they leave undetermined a dimcnsionlcss parameter (co s , in the 
notation of the first group) which was shown in p5[ to be crucial for determining the dynamics of the last few orbits, 
and (ii) they exhibit certain mathematical inconsistencies. By "mathematical inconsistency" , we mean here that 
both regularization methods defeat the very purpose for which they were used, i.e. to uniquely define a perturbative 
solution of Einstein's theory which depends only on two masses mi, rTT^- We recall that Einstein's field equations 
E^y = R^y — \R g^y = 8tt G T^ v exhibit the consistency feature that they imply the equations of motion of the matter 
{VT^y = 0) as a consequence of the (contracted) Bianchi identity VE^v = 0. To preserve this delicate consistency 
property when perturbatively solving Einstein's equations in some gauge it is necessary to use a regularization method 
which respects the basic properties of the algebraic and differential calculus of ordinary functions: such as the 
associativity, commutativity and distributivity of point-wise addition and multiplication, Leibniz's rule, Schwarz's 
rule (d^dyf = d u d^f), integration by parts, etc. The regularization methods used up to now at the 3PN level violate 
at least one of the basic properties of the standard calculus (for instance, Leibniz's rule is violated both by the 
regularization method of [p|J^] and that of 1 13 IJ]). 
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We know only one regularization method which formally respects the basic properties of standard calculus: dimen- 
sional regularization [|1 6| pL 9|| . It was invented as a way to preserve the gauge-symmetry of perturbative quantum gauge 
theories. Though our set up is classical, our basic consistency problem is to respect the gauge-symmetry (underlying 
the link between Bianchi identities and equations of motion) of perturbative general relativity. Dimensional regular- 
ization seems to be the ideal tool for generating a consistent perturbative solution of Einstein's theory depending only 
on mi, m2. In this letter, we shall show how to use dimensional regularization at the 3PN level in the ADM framework 
. Our three main results will be: (i) all the pole parts oc (d— 3) _1 (where d denotes the continued space dimension) 
appearing in intermediate contributions cancel in the total 3PN Hamiltonian iJ, (ii) the contributions to H which 
are quadratic in momenta are unambiguously determined by dimensional regularization to take the unique value 
(u>k = 41/24) which preserves global Poincare invariance J§|||, and (iii) the momentum-independent contributions to 
H are unambiguously determined and yield the following unique value for the heretofore undetermined dimensionless 
parameter entering 3PN dynamics: 

LJ s =0. (1.1) 



II. DIMENSIONAL CONTINUATION OF THE 3PN ADM HAMILTONIAN 

Let D = d + 1 denote the (analytically continued) space-time dimension. The ADM approach |ll| uses a d + 1 
split of the coupled gravity-matter dynamics and works with the canonical pairs {x a ,p a i) and (gij,Tt tJ ) (hj>k,... 
denote spatial indices taking (formally) d values; a = 1, . . . , N labels the particles). The dimensionally continued 
hamiltonian and momentum constraints read (in units where IQ-kGd = 1) 

V9 R = -jjj ('.'/<< 9H ^ - J^J i9ij 7T ij ) 2 ) + ^(ml+flPaiPaj)* S a , (2.1a) 

-2D 3 ^ =Y,9a°Pa 3 5 a . (2.1b) 

a 

Here 8 a = <5(x — x a ) (with J d d x8(x.) — 1), g l J = g tJ (x a ) (which will be seen to be perturbatively unambiguously 
defined and finite), and Dj denotes the c£-dimensional covariant derivative (acting on a tensor density). 

The application of the ADM formalism to the gravitational dynamics of N particles consists of four steps. Step (i) 
consists in fixing the gauge by demanding that gij and tt 1j have the forms (ADMTT gauge) 

9ij = A(d>) % + hj? , 7T« = 7T« (V k ) + 7T? T , (2.2) 

where 



d-2 



4/(d-2) 



A ^ = { 1 + W^T)V ' (2 ' 3a) 

n ij (V k ) = diV j + djV 1 - ^ ^ d k V k , (2.3b) 
and where the TT pieces hj? , ir^ T are transverse and traceless, i.e. satisfy dj f 7j T = = 8* 3 /? T with / = h or it 



Step (ii) consists in solving, by a perturbative expansion, the constraints (2.1 ) by expressing the "longitudinal" 



variables and tt v (i.e. V 1 ) in terms of the dynamical degrees of freedom (a;* ,p a i] hjj , 7r^ T ). This is done by means 
of a post-Newtonian (PN) expansion: — 0(2) + 0(4) + 0(6) + 0(8) + 0(io) > V* — Vfa + Vfa + Vfo . Here, the numbers 
within parentheses denote the formal orders in the inverse velocity of light. E.g., 0(2) ~ Gm / (c 2 r d ~ 2 ) . [Apart from 
this formal use of 1/c we set everywhere c = 1.] For instance, we have 

A 0(2) = -J2m a 5 a , (2.4a) 

a 

A ^) = -§E^+i^i)^)E m ^' ( 2 - 4b ) 

a a ^ ' a 

nr\ -AV.i _i_ I 1 —\ Ft. A/3 — _I . 



dj ^(1) (V) = A Vfo + I 1 - 2 ) d v % =-^PaiS a . (2.4c) 
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The solution of Eq. (2.4c) can be written as 

^ ' a a 

Step (iii) consists in expressing the total energy of the system E [<fr] (given by the surface integral at spatial infinity 
of —di 4>) in terms of the (particle and gravitational) dynamical degrees of freedom. This yields an Hamiltonian 

x aiPai, , ^tt) which describes the coupled dynamics of matter and gravity. The fourth and final step used in 
the ADM approach to the problem of motion of N point masses 21 23[f| consists in: (iv) eliminating the gravitational 



variables {hJ 3 T , ir^ T ) by perturbatively solving their field equations (as obtained from varying the Hamiltonian derived 

in the previous step) . This is again done by replacing hj^ and 7r^ T by the PN-expanded solutions of their field equa- 
tions. [We arc interested here in the conservative dynamics and use a time-symmetric (half-retarded half-advanced) 
Green function. See for the discussion of radiation-reaction effects in the ADM approach, and |23|]6],|7jl for the 

subtlety of the reduction of the higher time-derivatives appearing when solving (/i^ T ,7r^ T ) in terms of the matter.] 

These PN expansions start as: hj^ = + tiESy H , 7r^ T = vr^ TT + 7r^.N TT + ■ ■ ■ . Actually, the fact that the 

field equations for hj^ , 7r^ T must be satisfied implies that it is sufficient (at 3PN) to replace (hj?, 7r^ T ) by their 
leading contributions (hj^j, 7T(g\ TT )' These quantities satisfy the equations 

/ \ TT 

PaiPaj x d — 2 



^« = (-^^-2^*™.*u| - (^a) 



The superscripts TT in Eqs. (2.6) denote the application to a second rank tensor of the d-dimensional (spatially 
nonlocal) TT-projection operator. 

The matter Hamiltonian obtained after elimination of the gravitational variables can be simplified by using many 
integration by parts allowing one to successively eliminate 0(g), </>( 6 ), VL\ and Vfe. by using the elliptic equations 
they must satisfy. [Some details of this elimination process will be discussed in |2~0[| .1 For convenience, we have 
performed this elimination procedure in a form closely parallel to the one performed (when d = 3) in [|7j. Finally, the 
dimensional continuation of the matter-reduced dimensionally continued 3PN Hamiltonian reads (we do not write the 
d-dimensional versions of the non-problematic 1PN and 2PN Hamiltonians) 

#3PN = --4rV^+ [ d d x(h 1+ h 2 + h 3 ), (2.7) 

where 

( {A-d)[d-2f Q 2 (4-d)(d-2) c \^ 
kl = { 64(d-l)3 S ^ - 16(d-l) 2 S ^))\ m ^ 

d(d + 2)(3d-4) a3 W + 2KM-4) (d-4)(d + 2) \ p * A 

192(d-l)3 ^ + 16(d-l) 2 + 32(d-l) 2 S ^ ) \ ^-J" 

d(d + 6) a2 , d A 4-d \x-(pD\ d + A 



64(d-l) 2 "^ 8(d-l)" w 32(d-l) m3 " 32(d-l 



4 d „ . 2 3d 4 y 2 



+ 4(d3T)^ (7r (3) ) "2(M^^)) 
(4 — d)(d — 2) i tt (d- 2)(2d- 3) TT 1 TT x - p^ ni p a j 

4(d-l) 2 < P( 2 ).''V)J / V) l j 2 (d- l) 2 n 2 ).''^ 4 ) J""(4)y + 4 ft (4)y 2^ ^3 

- ((^)^3)) TT ) 2 - ^(^( 2 )-I)) TT ^ ~ ~ (^4) 2 • ( 2 - 8a ) 



5a 
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FIG. 1. "Three-loop diagram" representing the contribution (of order G 4 m? mi) proportional to J d d x 4>i SfcftSolij ^k^?40)y 
in the Hamiltonian, where hf^ij = —2c(d)A~ 1 (di4>idj(f)2) TT - The solid lines represent the "propagator" A -1 of 0, the 
helicoidal lines represent the "propagator" (A _1 ) TT of hj^p and the circles labelled 1 and 2 denote the sources mi 8i and 
mi 62- 
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■ TT 



„ fc 1 tt (d— 2)(3d— 4) TT 
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d + 2 



4(d-l) 



(2) /l (4)ij 



,fc 



(4)ij,k 



2(d- 1) 



(2)^(4)ij 



u TT 

J (4)ij,k ■ 



(2.8c) 



The 3PN Hamiltonian (2.7) is expressed in terms of the following quantities: 0( 2 ) defined by Eq. (2.4a), 0(4) defined 
by Eq. ( pb| ), V ( l 3) defined by Eq. @, tt^ = 7r <J '(V ( ^) defined by Eq. ( jOg ), h^. defined by Eq. ( |2~6a| ), as well as 



S'(4) and >S(4)jj denned by the following equations 



v / £ ^ ^ 



(4)y 



PaiPaj 



5 a • 



(2.9) 



Using the natural (propagator)* (source) structure of each building block of the Hamiltonian H, we can represent 
H by a sum of diagrams. For instance the last term in Eq. (2.8c) gives rise (after expanding it in powers of mi and 
TO2) to several diagrams, one of which (oc m\ m|) is represented (for illustration) in Fig. [I]. Note that this diagram 
contains three (classical) loops. More generally the most nonlinear (momentum-independent) contributions to ff„pN 
contain n loops. 



III. FINITENESS OF H A DM(d) AS d -> 3 

The perturbative consistency of dimensional regularization | fl9[ ensures that the formal d-dimensional PN-expanded 
solution g^uix^) given in Section [n| will directly generate, when d is analytically continued up to d — 3, a 3-dimensional 
solution of Einstein's theory if it is finite as d — > 3, i.e. if no poles proportional to (d — 3) _1 arise when d — > 3. One can 
analytically find out all the potentially "dangerous" terms (those which might give rise to a pole) in £i3PN, Eq. Q2.7|), 
by the following analysis. First, the d-dimensional solution of the basic equation Au a = —6 a is u a = — A _1 <5 a = k r„ , 
where k = V (^p) and r a = |x — x a |. [Though k depends on d we shall factor the various powers of k entering 

each contribution to i?3PN without ever expanding its dependence on d as d — > 3.] This allows one to write down 
explicitly 0(2), 0(4), an d The explicit expression of VA\ is then obtained by iterating the Poisson operator 

A" 1 using A- 1 ^ = ((A + 2)(A + d)). 



These results suffice to write down the part of hj^j, Eq. (2.6a), whose source is — m a 1 p a iPaj 5 a - We denote this 
part as hj^ ai - (with a = 1,2). The part of which does not depend on momenta will be denoted by hj^^-. One 
can control its crucial "ultra violet" (UV) behaviour as r a — > (which suffices to analyze the potentially dangerous 
terms) in the following way. First, let c(d) = (d— 2)/(2(d— I )) so that the source term for hj^ t - is the TT projection of 
— c(d) crij (0(2) , 0(2)) where <Jij (0, ip) = di<p djip. Inserting the solution 0( 2 ) = mi u\ + m 2 u 2 into <xy , and applying the 
TT operator (which projects to zero aij(ui, and cry (112, 112)), yields hj^.^ = — 2 c(d) mi m-2 A -1 (cry (tti, u 2 )) TT . 
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Operating by parts under the TT operator further yields hj^^ = +2 c(d) mi mi A _1 (wi <9yit2) TT . It is then easily 
proven that, near xi, /ij-^y is the sum of a regular (C°°) function and of a singular expansion in powers of r% 



obtained by inserting in the source term tx {u\ <9yit2) TT the Taylor expansion of <9y U2(x) near x = xi, and applying 



the operator A 1 to each term of the expansion. 

Using this local UV analysis we studied the convergence, near, say, xi, o f the integral giving Habm- The conclusion 



is that poles oc (d — 3) _1 can arise only from the third integrand /13, Eq. ( 2.8c ). More precisely, one finds that, when 
expanding each field entering /13 in powers of the masses m\ and 7712, and then expanding each partial integrand in 
powers of n there arise ten (and only ten) separate terms Ta, A — 1, . . . , 10, giving rise to poles as d — > 3. These 
pole-generating partial integrands have all the structure (for the poles generated by the integration near n = 0) 

T A = fc 4 miTO 2 rf~ 3d (c A iD(pi,pi) +c A2 (tii •pi)D(n 1 ,p 1 ) 

+ c A3 (m ■ pi) 2 D(ni,m) + caap\ D(ni,nifj , (3.1) 

where Cax, ■ ■ ■ , Caa are (i-dependent coefficients and where Dy = [<9ijr2 _d ] x _ , D(p, q) = p l qi , n\ = r^ 1 (x l ~ x\), 

and ri\ ■ p\ = n\ pu = Sij n\ p\ . These ten terms arise from the terms in the mass-expanded version of /13 with the 
respective structures: fa fa {dVi) 2 , dVi dVi hj^, Vi dV\ dhj^, fa dfa dfa hj^ v fa dfa dfa hj^ v fa dhj^ 2)1 dhj^, 

OSi dfa /i(4Q), OSi fa dhj^, dfa dhj^^ hj^, and dfa /i^i dhj^. Here we have suppressed, for readibility, the PN 
indices (2) (on fa, (3) (on V) and (4) (on S), but we left the indices (40) and (42) on /i TT . The explicit labels 1 or 2 refer 
to their sources oc Si or 62- Each term Ta c arries its full coefficient, as arising from the successive (mass-expanded 
and Taylor-expanded) terms in /i 3 , Eq. ( |2.8c ). 



One then computes the "local" contribution (in a ball around xi of radius £1, with < li <C 7*12 = |xi — X2I) of each 
"dangerous" term Ta- After averaging (in a d-dimensional sense) over angles, and integrating over the radius r\ , we 



find that each "dangerous" term locally yields, when d — > 3, the following "pole" contribution to the Hamiltonian (2.7) 



HY{d) = -\ n d k i mim 2 D(pi,pi)£f {d - 3) ^ , (3.2) 

where Qd is the area of the unit sphere in d dimensions, and where we defined 

C A{d)=CA 1 {d) + — r + w ^ y (3.3) 

We have explicitly computed the 10 x 4 = 40 ^-dependent coefficients CAi(d), . . . , CAi{d) as well as the 10 angular- 
averaged coefficients CA{d). Note that the value of each of these coefficients delicately depends on having correctly 
derived the d-dependent coefficients appearing in the Hamiltonian, as well as all the d-dependent factors arising from 
the A -1 and TT operators and from the angular averaging. The results are listed in Table I, each coefficient being 
expanded in powers of s = d — 3, e.g. CA{d) = ca(3) + sc' A (3) + 0(e 2 ). A first important result from Table I is that, 
as one checks, the sum of the ten 3-dimensional averaged coefficients 0,4(3) exactly vanishes: 
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5>43) = 0. (3.4) 



A=l 



In view of Eq. (3^2) this proves that the total pole part, as d — > 3, of H{d) vanishes. In other words, dimensional 
continuation up to d — 3 gives a unique, finite value for H^f^. [This result confirms the finding of that all the 
"Riesz poles" appearing in the Riesz-based implementation of Hadamard's partie finie finally cancel. However, this 
finding was not sufficient for being able to unambiguously compute the "regularized" value of -£^3^^ •] 

IV. DIMENSIONAL-CONTINUATION DETERMINATION OF THE HAMILTONIAN 

Having shown that H^^ A (d) has a unique, finite value as d — ► 3, it remains the crucial task of explicitly computing 
this finite val ue. L et us first focus on the /13-generated part of the Hamiltonian, i.e. on H^(d) = J d d x /13(d), with /13(d) 
given by Eq. ( 2.8cj ), and on the finite limit: lim,j— >3 H^{d). Let us compute the difference AH3 = lim,j_*3 H^{d) — i?P JS , 



where i?.p js = RH \J d 3 xh§ JS ] is the Riesz-implemented Hadamard regularization (RH), as defined in HJl^M, of the 
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TABLE I. Expansion in powers of e = d — 3 of the d-dependent coefficients CAi(d), . . . , CA4,{d) and CA{d) entering the fO 
"dangerous" terms Ta, A = 1, . . . , 10. 



.4 



CAl 



CA2 



('A?, 



CA4 



(-A 



315 3711 315 267 21 1013 

1 e 1 e 1 e 

1024 4096 2048 1024 512 10240 

_L3___263_ 5L_iL ii. .5! JL 2^. 7 1649 

' 256 ~ 1536 6 ~ 128 ~ 256 6 512 + 256 £ 512 + 3072 6 ~ 128 ~ 7680 6 

7 487 37 91 29 665 49 49 7 601 

+ — - — £ — TT^T — TTT-T £ ~ - — £ — — — - -f" T^TT £ -— — — + 



128 6144 256 384 1024 6144 1024 768 2560 153600 

15 191 5 57 1 53 



512 2048 512 2048 256 5120 

__5___37_ 5 57_ _5_ _37_ 1 33 

" 64 ~ 256 6 ~ 128 ~ 512 £ 128 + 512 6 ~ 32 ~ 640 6 

43 7 31 221 73 61 17 263 49 5467 

1536 + 256 £ 768 + 1536 £ ~3072 ~6144 e _ 3072 _ 6144 £ 1280 + 76800 6 

65 421 25 35 _ _5_ _ _37_ _^L:JtL _L 15 

768 + 1536 6 ~ 384 ~ 256 £ ~ 1536 ~ 3072 6 1536 + 1024 6 16 + 64 6 

__65_ _ _161_ _25_ _5_ _5_ _L7_ _JiL _LL _i L 

~768 ~T536 e 384 + 768 £ 1536 + 3072 £ ~ 1536 + 3072 £ ~16~64 £ 

13 23 19 49 7 11 23 257 21 5423 

■ £ — — — — — £ — — — + ■ 



3072 1536 1536 768 6144 4096 6144 12288 2560 153600 

13 5 19 5 7 61 23 55 21 2903 



3072 768 1536 128 6144 12288 6144 4096 2560 153600 



3-dimensional integral of the integrand h^ 3S given by Eq. (A9c) of Note that each term in /13(d), Eq. (2.8c), tends, 



when d — > 3, to a corresponding term (with the same coefficient in the limit) in h^ is . The comparison between the 
two integrals H3 (d) and i?.p js is best done by: (i) separating the full space integral into two local integrals (restricted 
to balls of radii l\ and £ 2 around xi and X2) and a global integral over the rest of space, and (ii) decomposing each 
term of each integrand into a regular part (which is absolutely convergent in d — 3) and a singular part (which is not 
absolutely convergent in d = 3) made of a finite number of terms of the form (say, near Xi) r^ Ci 1 ...i k ■ ■ ■ n\ k . It 
is then seen that the difference AH3 is simply given by the difference between lim^a H 3 oc smg (d) and RH [H^^^], 
involving only the local (rj < i\ or r2 < £2) integrals of the singular parts of the integrands. Next, one checks that all 
local singularities which do not correspond to poles in d — 3 are given, in the limit d — > 3, the same regularized values 
in both regularization methods. However, the situation is different for the local ("logarithmic") singularities oc r^ 3 
in d = 3, which give rise to poles in d 7^ 3. The complete list of these "dangerous" local singularities has been given 
above: they are the ten terms Ta listed in Table I. [Note that they are all quadratic in p\. Here, we focus on the UV 
behaviour near Xi; one should add the similar terms obtained by exchanging 1 «-* 2.] After averaging over angles, the 
RH regularization of the d — 3 limit of each local singular integral is oc ca(3) ln(£i/si) where si is the scale entering 
the definition of the RH operation. After summing over A = 1, . . . , 10, the result (O) shows that RH (H^g*) = 0. 



On the other hand, dimensional continuation gives for the value of H Mng (d) the sum over A of the values (t 



Finally, we conclude, using ca(3 + e) = c^(3) + e c' A (3) + C(e 2 ), Eq. (3.4), and the 3-dimensional limits = 47r, 
k(d = 3) = 1/(4tt), that (r 12 = |x x - x 2 |) 

AH 3 = Bmff,(d) - i/ 3 DJS =-\(fl c ^ 3 )) PUPV 9 ^ + 1 " 2 ■ 

We have similarly evaluated AH t = lim d ^ 3 Hi(d) — ffj DJS , for i = 1,2, where ifj(d) = / d d xhi(d) and ijp js = 
RH [/ d 3 x /if JS ] . Here the formal expressions of h\(d) and /i 2 (d) do not seem to correspond to the expressions of hY 3S , 



G 



given in Eqs. (A9a), (A9b) of [Q. However, we can use the properties of dimensional regularization to transform 
h\{d) and h 2 (d). We have checked that, by using suitable integration by parts (which are allowed in dimensional 
continuation), one could transform hi{d), i = 1,2, into new expressions h'^d) which are term- by-term similar to hf^, 
in the same sense that /13 (d) was similar to JS . An analysis (using the tools given above) of the UV behaviour of 
h'^d) then shows that they contain no "dangerous" terms 0(r 1 x ^) with A(3) = 3. The reasoning of the previous 
section then shows that AH\ = AH2 = 0. 

Finally, we conclude t hat the finite value of the total Hamiltonian is given by: lim^s H^^ 1 (d) = i/fp^ + AH3 
with AH3 given by Eq. ( |4.l| ). The two main consequences of this result concern the two dimensionless parameters 
that were left undetermined by the original Riesz-Hadamard regularization of the 3PN Hamiltonian [|| : uj^ and uj s . 

First, from the definition (see, e.g., Eq. (13) of j^[) of the dimensionless "kinetic" parameter u>k and the fact that 
i?.p js corresponds, by convention, to uj^ = 0, the result ( f4.l| ) is equivalent to saying that dimensional continuation 
uniquely determines the value of u>k to be (using (47r) _1 = AG) 



dim rcg 



64^4(3). (4.2) 



A=l 
41 



The sum of the c' A (3) read from Table I is checked to be ^||g . The dimensional continuation prediction (L2) therefore 

yields w^ lmics = As shown in J§| (see also p|JlO||) this value is the unique, correct value which ensures that the 
regularized Hamiltonian admits a nonlinear realization of global Poincare invariance. We find it truly remarkable that 
the use of dimensional regularization in the ADM formalism (which violates from the start manifest Poincare invariance 
by splitting space and time and by using a non Poincare-invariant gauge) uniquely predicts, after a long calculation 
involving many intermediate complicated d-dependent coefficients and finally summing 50 different contributions of 
Table I (indeed c' A (3) is a combination of c' Al , c' A2 , c' A3 , CA2, and ca3 for A = 1, . . . , 10), the unique, physically 
correct value of ojk- We interpret this as a very strong confirmation of the mathematical (and physical) consistency 
of dimensional regularization, and, in particular, of its property of perturbatively respecting gauge symmetry (i.e., in 
our case, diffeomorphism invariance). 

Second, as AH3 affects only the part of the Hamiltonian which is quadratic in the momenta, we conclude that 
the heretofore undetermined "static parameter" lu s is simply: lu s — (as it was in the "reference" Hamiltonian 
jj" DJS ). Note that this simple conclusion was obtained after a detailed analysis of the UV behaviour of all the terms 



(2.7). In this analysis, the ADMTT gauge played a very useful role in suppressing many of the stronger UV 



divergencies that occur in a harmonic gauge. For instance the contribution I\{d) — J d d xcf>i dkhj^^ dkhj^^ where 
hj^,^ = — c(d) A~ 1 (di4>djcj)) TT , which corresponds to the diagram depicted in Fig. [l], turns out to be completely "non 

dangerous" (and even locally integrable), while its harmonic-gauge analog would contain Uij = —c(d)A~ 1 (di4>dj(j>) 
instead of an d would give rise to several (gauge-spurious) "dangerous" integrals generating pole parts (and 

logarithms in d = 3). [The existence of poles at 3PN, in harmonic coordinates, was first mentioned in [|| and 
has been confirmed by the explicit results of fjcj] .] Let us also note that dimensional regularization leads to a 
very simple and consistent treatment of "contact terms" . Starting from the best available definition of dimensional 
regularization, which is in the Fourier domain (e.g. (l^]), we have checked that the value of, e.g., hid) — J d d x </> 4 (x) <5i 
is </> 4 (xi) = (m 2 U2(x!)) 4 . This simple "threading" property of contact terms extends to all the perturbative contact 
terms arising when formally solving Einstein's equations. Here, dimensional regularization differs from some variants 
of Hadamard's partie finie regularization in which Pf {4> ) 7^ (Pf (</>)) . 

To conclude this letter let us mention some of the physical applications of our results. The most important is that 
our unique determination of the last missing parameter in the 3PN Hamiltonian allows one to make full use of the 
3PN accuracy and, in particular, to analytically estimate the characteristics of the last stable orbits (LSO) of binary 
black holes. For instance, using the best available PN-resummation method, namely the (Pade-resummed) effective 



one-body approach [£5],[l5 26 1, one finds that the LSO binding energy of two non-spinning equal-mass black holes is 
1.67%(mi + m2)c 2 . More generally, our work allows one to extend to 3PN the approach proposed in [ p5| , namely to 
start the numerical simulation of the coalescence of two black holes (without or with spin |2q |) just after they cross 
their LSO, so that they are less than one orbit away from coalescence. 

Our work contains another important lesson for numerical relativity: our dimensional-continuation calculation of 
i?3 and u>k shows (as a corollary) that the "conformally flat" truncation (i.e. setting hj^ to zero) which is of general 
use in numerical relativity, is "inconsistent" in that it ruins the delicate cancellation of poles and leads to a formally 
infinite 3PN Hamiltonian. [In addition, it violates global Poincare invariance.] As the physical meaning of the 
3PN pole cancellation is that the perturbative algorithm unambiguously determines the unique, physically correct 
solution of Einstein's equations depending on m\ and mi which matches two black holes |20| ], the lack of uniqueness 
associated to the arbitrariness in subtracting the 3PN pole generated by the conformally flat truncation probably 
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explains why current initial data simulations [|27],|28| significantly differ from the resummed analytical estimates of 
LSO characteristics. 

Finally, our work having convincingly established the perturbative consistency of dimensional regularization, we 
strongly recommend to use it to uniquely determine the 3PN contribution to the gravitational wave flux, whose 
determination might otherwise be marred by spurious ambiguities. Let us also mention that it would be nice to 
independently confirm our results by two different calculations: (i) to use dimensional continuation, instead of a 
variant of Hadamard regularization Jl3| , ^4[ , to compute the 3PN dynamics in harmonic coordinates, and (ii) to use 
dimensional continuation to recompute uj s by studying the motion of a test mass around a Schwarzschild black hole 
to the first "post-test-mass" approximation (0(fj,/M) corrections). 
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